Entanglement and dynamics of spin-chains in periodically-pulsed magnetic fields: 

accelerator modes 
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We study the dynamics of a single excitation in a Heisenberg spin-chain subjected to a sequence 
of periodic pulses from an external, parabolic, magnetic field. We show that, for experimentally 
reasonable parameters, a pair of counter-propagating coherent states are ejected from the centre of 
the chain. We find an illuminating correspondence with the quantum time evolution of the well- 
known paradigm of quantum chaos, the Quantum Kicked Rotor (QKR). From this we can analyse 
the entanglement production and interpret the ejected coherent states as a manifestation of so-called 
'accelerator modes' of a classically chaotic system. 
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There is considerable interest in the fidelity of quan- 
tum state transmission and entanglement measures in 
spin chains because of their relevance to quantum in- 
formation applications. In Q state transmission in a 
Heisenberg chain was investigated. In |2| it was shown 
that such a chain, in the presence of an external, static, 
parabolic magnetic field, can give perfect transmission of 
coherent spin states of appropriate width. Obtaining co- 
herent states of specified widths, represents a technical 
challenge, though. 

Here, we investigate the dynamics of a Heisenberg spin- 
chain subjected to short, time-periodic pulses from an 
external parabolic magnetic field. We find that this pro- 
vides an effective technique for generating well-defined 
coherent states, starting from a single excitation at the 
centre of the spin chain. The key to the analysis is that 
we note, for the first time, a close correspondence be- 
tween the time-evolution of a Heisenberg chain and that 
of the well-known chaotic system the Quantum Kicked 
Rotor (QKR) 0, Q in its quantum resonance regime. 
Our additional parabolic external field extends the cor- 
respondence between the Heisenberg spin-chain to the 
non-resonant QKR. The non-resonant, chaotic, QKR has 
been well-investigated experimentally with cold atoms in 
optical lattices Q ■ The QKR in the resonant regime has 
also been investigated experimentally Q. 

There is also much current interest in the interface 
between quantum chaos and quantum information 
In some studies of entanglement measures, the quantum 
chaos is generated by extrinsic disorder , in others with 
a clean but chaotic Hamiltonian The question of 

whether chaos aids or hinders entanglement generation 
relevant to quantum information applications has not yet 
yielded a clear answer ^IQj. In it was shown that a 
class of kicked Ising-type chains have quantum behaviour 
related to those of one-body 'image' systems with a well- 
defined classical limit, which can be chaotic or integrable. 

But, to our knowledge, the correspondence between 
the dynamics of the QKR, a leading paradigm of quan- 
tum chaos and the Heisenberg chain, a system of such 
key interest in quantum information, has not been noted 



or exploited previously. We show it means that with 
the pulsed parabolic field, we can employ certain 'text- 
book' 01 expressions found for the QKR and the Stan- 
dard Map to describe the entanglement properties. It 
means also that we see not only generic forms of quan- 
tum chaotic behaviour in the spin-chain like exponen- 
tial localization (analogous to Anderson localization seen 
in disordered metals) but will also generate phenomena 
(such as 'accelerator modes') due to additional correla- 
tions specific to a 'clean' chaotic system -and to the QKR, 
in particular. We show that, as we can remain in the one- 
excitation sector with this Hamiltonian, there are new 
possibilities for quantum information applications from 
the entanglement properties of the accelerator modes. 
We consider a time-periodic Hamiltonian of the form: 
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where Hhc = -i E„ ^r" ' ^^"^^ ~ E„ ^cr? is the Hamil- 
tonian for the Heisenberg chain studied in 0. Tq is the 
period of the pulses; Bq is the amplitude of the applied 
parabolic magnetic field; the length of the chain, TV > 100 
in the present work. The spin-transmission properties for 
the time-independent part were investigated in I'l : for a 
non-zero static field, where hB 3> kT one may restrict 
the study to the single excitation regime (ie restricted to 
the basis of states |s), which have a spin-up at a single 
site s on the chain but all other spins down (along —z). 

In [2 it was shown that the eigenstates of Hhc , 
are delocalized along the chain ie 
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Using the eigenstates, an ana- 



lytical form for the time-evolution operator Uhc{t, 0) = 
eyi?p{—j^Hhct} may be obtained in the single-excitation 
basis: 
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where we disregard the uninteresting overall phase due 
to the uniform static field B (or formally set 2BTo = 
27r). For the periodically-pulsed system described by the 
Hamiltonian in , the matrix elements of the one-period 
time-evolution operator U{Tq, 0), in the single excitation 
basis, are: 
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In FigQ] we show the effect of repeated application of 
© and Q on the state ip{t = 0) = |rio) (a state initial- 
ized on a site at, or very near, the centre of the chain. The 
spin-amplitude spreads out into a very irregular 'chaotic' 
distribution around the site Uq. But most strikingly, we 
can see a pair of counter-propagating spikes, 'hopping' 
around 2'k/Bq ~ 94 spin-sites each consecutive period. 
In order to analyse this behaviour, we re-examine the 
form of [//"/(TcO). 

For large TV, it is easy to see that becomes the dis- 
cretized version of an integral: as — > oo, (r|[/'"^|s) — > 
Fra where 



Frs — / [cos(r H- s — l)x + cos(r — s)x] (5) 
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The X — ■n{m — 1)/N is nominally a position coordi- 
nate but in fact represents motion through the subspace 
of I to), the eigenstates of Hhc- Eq. © can be compared 
with the time-evolution operator of one of the most ex- 
tensively studied system in quantum chaos, the Quantum 
Kicked Rotor (QKR) 0. 

The QKR corresponds to the Hamiltonian H = ^ 

K cosxJ^nH^ ~ '^)- where K is the kick strength and 
T = 1 is the kicking period. The equations of motion for 
its classical limit produce the 'Standard Map', the text- 
book paradigm of classical Hamiltonian chaos T^]. The 
time evolution operator of the QKR, {/Q^^ is generally 
given in a plane wave basis jj]; ie 

(n\UQ^^\l) = e-''''^/^ i''-^J^_i{K/K) where J denotes 
an ordinary Bessel function. Note that in the so-called 
quantum-resonance regime, where h = 27r, (also stud- 
ied in cold atom experiments |Q|) we obtain simply 
{n\U'^^^\l) = i"'-^Jn-i{K/h). In the experiments, h 
represents an effective value obtained from the optical 
latticeparameters and kick period T, typically of order 

However, since parity is conserved, the basis states of 
the QKR are generally symmetrised, ijjf^ = l/\/2[\l) ± | — 




FIG. 1: Time evolution of the state ^p{t = 0) = |no) (ie initial- 
ized with a single spin-excitation at the centre of the chain)- 
showing the effect of accelerator modes. P{s) represents the 
probability of finding the excitation at site s. The accelerator 
modes are the 'spikes' at the leading edge of the distribu- 
tion. They correspond to a counter-propagating pair of co- 
herent states ejected from the centre. We take Bq = 1/15 and 
2JTo = 100. The upper line is at f = 3To; the lower curves 
correspond to consecutive periods jTo with j = 4,5,6... as 
numbered. The dotted line indicates the form of the parabolic 
field (scaled by a constant factor) which is pulsed on/off every 
period at times t = JTq. The accelerator modes represent over 
25% of the total probability; they advance an equal distance 
(shown below to be 2-n /Bq ~ 94 spin sites) each period, and 
after just 3 pulses are well separated from the central, 'chaotic' 
remnant. 



I)] with I — 0, 1,2, ...M, where M is where we truncate 
our basis, and |/) = l/\/27rexp{iZa;}. 

It is evident that, to within an overall phase, the form 
of the integral in is equivalent in form to the 'kick' 
part of the QKR operator. 
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if 2JTq = K = I3h\ above, \r— 1/2) indicates plane waves 
shifted 1 
Further, 



shifted by half a quantum ie \r — 1/2) = — i=e*^'' -^/^^'^ 



\i^U/2) ^ ^"-'Jr-siP) (7) 



provided we neglect terms which are only significant at 
the edges of the basis (ie r w 1 or r « N) . Note that for a 
spin-chain on a ring, these 'edge' corrections are entirely 
absent and the time evolution of a Heisenberg chain ring 
is entirely equivalent to the kick part of the QKR. The 
addition of the parabolic-field 'kick' completes the anal- 
ogy with the full QKR: the fact that the parabolic field 
is now the S — kick term, while for the QKR the kinetic 
energy provides the time-independent term, represents 
only a simple re-scaling of the parameters in U{Tq, 0). 
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Note that all our numerics employ and not the 
QKR form in {Tj) above. Nevertheless, we found the spin- 
chain dynamics are sufficiently analogous that we can 
simply make the substitutions: (a) Bq = his an effective 
value of Planck's constant, (b) 2JToBq — K, where 
K is the stochasticity parameter which fully determines 
the classical dynamics of the kicked rotor; eg if > 4.5 
indicates fully chaotic dynamics. We can then directly 
apply the well-known QKR results to the spin dynamics: 

1) Short-time diffusion of spin- excitation: For K > 
4 the QKR initially displays diffusion of momentum, 
(p^) « D{K) t, which follows the classical behaviour: 
the form D{K) « K'^/2{1 - 2J2{K) + 2J2{Kf} [ll is 
obtained from a study of the classical diffusion for the 
Standard Map. The first term in the expression for D 
represents uncorrelated diffusion (equivalent to a random 
walk); the J2{K) terms result from short-time correla- 
tions present in the classical kicked rotor. We obtain 
analogous behaviour in our pulsed spin chain; substi- 
tuting 2JTqBq = Kg, and starting out with a spin ini- 
tialised on a site s — sq anywhere on the chain, wc find 
((s — so)^Bq) « D{Ks)t. The diffusion of spin excita- 
tion obtained from the time-evolution operator in Eq|31is 
compared, in Figl^^a), with the form of the classical diffu- 
sion rate of the QKR. For Kg = 10, the correction due to 
short-range correlations is large, but for Kg = 5, we find 
D K. Kl/2 since J2(5) « 0. After the so-called 'break- 
time', t* ^ (Ks/Bq)'^ the quantum diffusion no-longer 
follows the classical behaviour and the spin-excitation 
stops spreading onto neighbouring sites. 

2) Dynamical Localization: The Floquet states (eigen- 
states of U^^^{T, 0) ) of the QKR are known to be de- 
localized in the resonant regime; hence, as investigated 
in 0, a spin-excitation can spread (imperfectly) from 
one end of the chain to the other. In the non-resonant 
case, however, the corresponding Floquet states exhibit 
Dynamical Localization, a phenomenon analogous to An- 
derson Localization in a disordered metal (J]. Hence, for 
the spin-chain in pulsed parabolic field, a single excita- 
tion anywhere in the chain will spread up to a maxi- 
mum and will 'freeze' (localize) at the break-time, t* . 
For t > t*, 

Plods) --^exp{-2\s~so\/L}. (8) 

The quantum locahzation length L ~ (2JTo)^/4. Fig. 
|2Ib) demonstrates the typical probability distribution of 
a dynamically localized state, obtained by evolving a 
state initially in an arbitrary site. 

3) Entanglement: Our knowledge of the time evolution 
and spreading of the spin distribution enables us to es- 
timate some entanglement measures; in particular, the 
Q-measure and the Concurrence The former is 
a measure of the global entanglement for a pure (multi- 
partite) state: 

Q(i^))=^(i-|:Ki^)4(i-i) 



where the ak are the amplitudes when the state is pro- 
jected on to the basis of spin sites, {|fc)}. We note the 
relation of this formula to the well known Inverse Par- 
ticipation Ratio, R = for a localised system 
R~2L. 

We assign a measure to the bipartite entanglement of 
two separated sites i and j with the Concurrence. Given 
a pure state in the single excitation basis Cij = 4|Q;i||Q;j| 
|17j . With the exponentially localized form, ((SJ: 

^exp|-|-(|i-sol + b'-Sol)| 

if the separation between two sites dij = \i — so\ + \ j ~ 
Sol > 0, the maximum of the concurrence occurs ai L — 
2d, with the value ■ 

4) accelerator modes: For values oi K ^ 2n, transport- 
ing islands of stability re-appear in the classical phase 
space. Their phase-space area is ^ 1/10 and hence can 
support quantum states if h < 1/10. Their effect has 
been experimentally observed for cold cesium atoms in 
pulsed standing waves of light |0| . Because of the rela- 
tively large effective values of ft > 1 in the atomic exper- 
iments, their effect there was diluted by the chaotic 'sea' 
of trajectories; nevertheless, the accelerator modes (and 
Levy fiights due to their environs) manifested themselves 
as an enhancement of the diffusion for X ~ 27r. 

In the spin-chains there is no evident bar to a low 
h = Bq : it simply requires a weaker magnetic field. Ac- 
celerator modes, in fact, are stable for a broad parameter 
range K = a27T where a = 1.03 — 1.10 so there is no need 
for very precise fine tuning of the parameter 2 JTq. The 
associated islands of stability span a width (of spin sites) 
of order As ~ j^2n/h, (about 10 sites for Fig. ^ so 
there is no need to prepare the initial state exactly at the 
centre of the chain either, in order for the initial state to 
overlap strongly with the accelerator islands. 

We find useful coherent states are obtained provided 
Bq < 1/5. The width of these states is simply deter- 
mined by the effective h : we verified numerically that 
their form is given by ip{s,t = JTq) ~ Aj exp{— i?Q(s — 
Sj)^} where sj = 2Tr j/h and j is the pulse number. We 
find that, for j ~ 1 the coherent state pair represents 
about 30% of the total probability. For finite ft we expect 
the amplitude to gradually 'tunnel' out of the accelera- 
tor island: for ft = 1/10 we estimated the decay numeri- 
cally \Aj\'^ ~ |Aj=ip exp{— j724}, so even after 30 pulses 
there is a substantial amplitude. For ft = 1/15 (shown in 
Fig. ^ the island can in fact support more than one 
eigenstate, leading to a slight oscillation in its ampli- 
tude. It only requires 2 — 4 pulses to cleanly separate 
the travelling states from the chaotic remnant, though, 
so even ft ~ 1/5 would give coherent states with proba- 
bility > 20%. These could then be taken into a static, 
parabolic field region and transmitted onwards with per- 
fect fidelity (2]. 

For an actual realization we might suggest Bq ^ 1/10; 
in an experiment, Bq = BQdt where St <C Tq is the 
pulse duration. If our maximum magnetic field Bmax is 
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FIG. 2: (a) Shows that, the rate of spreading of a spin exci- 
tation initially at an arbitrary site |so) in the chain, at short 
times, is determined solely by the 'classical stochasticity pa- 
rameter' Ks = 2JToBq. At longer times, the spreading sat- 
urates. The straight line indicates the behaviour expected of 
the classical chaotic kicked rotor. The smaller the effective 
value of Planck's constant, Bq, the longer the spin spreads 
linearly at the 'classical' rate. Note that in the accelerator 
mode regime however, the increase would be quadratic in 
time, not linear, (b) Shows that at sufficiently long times 
the spin probability distribution saturates into an exponen- 
tially localized form. This is the analogue of Anderson Lo- 
calization in a disordered metal. P{s) ^ exp{— 2ls — sq\/L} 
so ln{P{s)) ~ 2/L takes the characteristic 'triangular form'. 

The straight lines correspond to L = — for Ks — 10 and 

agree well with the numerically calculated distributions. 



of order 0.1 Tesla, B^ax ~ [N/2YBq - 10"*^ au , for 
AT - 100 - 1000, implies Bq - 1Q~^° au so tlie pulse 
duration 5t ~ 10^ au ~ 25 ns. As we have made a split- 
operator approximation in Eq. Q), we require 2J5t is a 
small phase, whereas 2 JTq is significant. This constrains 
2J ~ 10^® ~ 1 — 10 MHz, so we are in a weak spin- 
coupling regime, for this choice of parameters. 

We now propose an apphcation of the accelerator 
modes in context of quantum communication. It has 
been noted that encoding quantum information in Gaus- 
sian wave-packets of excitation (where the presence and 
absence of the wave-packet depicts logical |1) and |0) 
states of a qubit) is a useful way of transmitting it down 
a Heisenberg coupled chain of spins '5, To"]. However, 
precisely how to create a superposition of the presence 
and absence of a Gaussian wave-packet remained un- 
clear. The work presented here suggests, in fact, that 
after 3-4 kicks one creates a superposition of two Gaus- 
sian wave-packets traveling in opposite directions and a 
exponentially localized state in the middle of a Heisen- 
berg spin chain. If one measures the exponentially lo- 
calized part and does not find the spin excitation there 
(which can have nearly 30% probability of occurrence for 
appropriately chosen parameters), then the spin chain is 
projected on to a superposition of oppositely traveling 
Gaussian wave-packets of excitation. If one denotes the 
left and right traveling wave-packets as \Gl) and \Gr) 
(which are each a Gaussian distribution of spin up), and 
|0l) and jOfl) as the absence such a wave-packet (which 
are each an all spin down state at the sites where \Gi,) 
and \Gr) would otherwise be) then the maximally entan- 
gled state \Gl)\^r) + \^l)\Gr) is effectively created. This 
state can now be distributed among well separated par- 
ties by switching on a constant parabolic field (instead 
of the kicks) in which a Gaussian wave-packet can travel 
for a significant distance without distortion . 
In conclusion: we have demonstrated that there is a 
close correspondence between a spin-chain in a pulsed 
quadratic field and the well-known chaotic kicked rotor; 
we find potential applications in quantum information 
processing. 
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